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ABSTRACT 

Vega and Fomalhaut, are similar in terms of mass, ages, and global debris disk properties; therefore, they 
are often referred as "debris disk twins". We present Spitzer 10-35 /Ltm spectroscopic data centered at both 
stars, and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The 
properties of the warm excess in Vega are further characterized with ancillary photometry in the mid-infrared 
and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many 
similar properties with the one found around Fomalhaut. The emission shortward of ^30 fim from both warm 
components is well described as a blackbody emission of ^170 K. Interestingly, two other systems, e Eri 
and HR 8799, also show such an unresolved warm dust using the same approach. These warm components 
may be analogous to the solar system's zodiacal dust cloud, but of far greater mass (fractional luminosity of 
~10~ 5 - 10~ 6 compared to 10~ 8 - 10~ 7 ). The dust temperature and tentative detections in the submillimeter 
suggest the warm excess arises from dust associated with a planetesimal ring located near the water-frost line 
and presumably created by processes occurring at similar locations in other debris systems as well. We also 
review the properties of the 2 pm hot excess around Vega and Fomalhaut, showing that the dust responsible 
for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may 
arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm 
and cold belt is rather large with an orbital ratio >10 in all four systems. In light of the current upper limits 
on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying 
planetary architecture, and suggest that multiple, low-mass planets likely reside between the two belts in Vega 
and Fomalhaut. 

Subject headings: circumstellar matter - infrared: stars, planetary systems - stars: individual (Vega, Fomalhaut) 
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1. INTRODUCTION 

Debris disks were discovered by IRAS as infrared excess 
emission arising from systems of particles analogous to en- 
hanced Kuiper belts. They are tenuous dusty disks sustained 
by cometary activity and planetesimal collisions, which initi- 
ate casca des of further collisions t hat break bodies down into 
fine dust (Dominik & Decin 2003). The dust is lost relatively 
quickly through photon pressure, Poynting-Robertson (P-R) 
or stellar-wind drag forces, generally in a time scale much 
shorter than 10 4 yr. Thus, a debris-disk-generated infrared 
excess requires the presence of colliding planetesimals, and 
a larger object (can be as small as Pluto-size) to stir them, 
i.e., some form of planetary system (Wvatt 2008). The large 
emitting area of debris makes these disks detectable through 
infrared/sub-millimeter thermal emission or optical scattered 
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light, providing insights into the nature of unseen parent-body 
populations and massive perturbers around other stars. 

The dust around a host star absorbs stellar radiation and 
re-emits in the infrared at equilibrium temperatures, balanced 
between absorption and emission. Consequently, the spec- 
tral energy distribution (SED) of the infrared excesses and 
its relation to the dust temperatures can be used to infer the 
number of emitting grains and their distances from the heat- 
ing star. For example, excesses in the near-infrared near 2 
/xm should be dominated by dust at ^1500 K, excesses in 
the mid-infrared near 24 /Ltm should be dominated by dust at 
^120-150 K, while excesses in the far-infrared are mostly 
from dust at ^50 K. Furthermore, excess emission at shorter 
wavelengths also contributes excess at longer wavelengths as 
Rayleigh-Jeans (a steep function of wavelengths). 

Identifying excess emission around a star requires a good 
knowledge of stellar photospheric properties and extrapola- 
tion to long wavelengths. Positive identifications of excesses 
at long wavelengths are easier than at shorter wavelengths 
where the host star dominates, unless the signal of the star 
can be filtered out through interferometric techniques. There 
is also a concern whether excesses at mid-infrared and far- 
infrared are tracing separate components since the majority of 
debris disks are unresolved. For nearby systems, modern fa- 
cilities like Spitzer and Herschel have sufficient resolution to 
resolve the detailed structures of the disk and reveal the com- 
plexity in disk structures in a few cases. The identification 
of a warm excess in these resolved systems requires precise 
subtraction of the stellar photosphere in the resolved images. 
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which has b een done for the Fom alhaut ( Stapelfel dt et al 
e Eri dBackman et al.ll2009l) . and HR 8799 dSu et al 



2004), 



2009) systems. Ironically, some of the nearby resolved sys- 



2 



tems are saturated in the mid-infrared, making the recognition 
of such a component very challenging. 

In this paper, we present mid-infrared spectroscopic studies 
centered at the two nearby A-type stars Vega and Fomalhaut 
obtained with the Spitzer Infrared Spectrograph (IRS) instru- 
ment and we identify warm excess emission in the vicinity of 
the star, clearly separated from the much brighter cold plan- 
etesimal belts that dominate the far-infrared emission. We 
suggest the presence of a planetesimal belt near the water- 
frost line of the Vega system for the first time. We com- 
pare its properties with the warm excess around Fomalhaut 
(Stapelfeld t et alJ [20041) . We discuss the properties of the 
warm excesses complemented with resolved infrared and sub- 
millimeter images of both sy stems. Two other sp atially re- 
solved debmsystems, e Eri (Backman et al. 20091) an d HR 
8799 (: Su et al.ll2009h . are also known to possess a similar 
warm belt that is fainter and clearly separated from the more 
prominent cold belt using a similar approach. These warm 
components may be analogous to our asteroid belt, but of far 
greater mass. We discuss the implications and origins of this 
two-belt configuration in light of the similarity in dust temper- 
ature distribution found around unresolv ed debris systems be - 
tween solar-like and early-type stars by Mora les et al.lfl2.Ql ll) . 

The paper is organized as follows. The observations and 
data reduction are described in Section 2, including both 
Spitzer IRS spectroscopy and Herschel PACS imaging. Pho- 
tospheric determination using ancillary data is given in Sec- 
tion 3.1 where we conclude that no significant excess is found 
from 2.2 to 8 pm for both systems. Using the PACS 70 and 
160 pm images, we estimate the excess flux of the unresolved 
source centered at the star position for Vega in Section 3.2 
and Section 3.3 for Fomalhaut. We construct the SEDs of the 
unresolved excess components along with other mid-infrared 
and submillimeter measurements, and suggest that they arise 
from a planetesimal belt placed near the water-frost line in 
Section 4. In Section 5, we discuss the implication of the 
two-belt systems on the underlying planetary configuration, 
and provide a possible explanation for the 2 /im excess in the 
Vega system. Conclusions are given in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer IRS Spectroscopy 

To avoid saturation, we only used data taken in the IRS 
high-resolution (R~600) modules (short-high (SH): 9.9-19.6 
pm and long-high (LH): 18.7-37.2 pm). The sizes of the slits 
are4."7xll"3 and ll"lx22."3 for the SH and LH modules, 
respectively; a significant fraction of the light from a point 
source is outside the slit. Standard slit loss correction for a 
point source was applied to the extracted spectra; therefore, 
the part of the spectrum where a point source dominates the 
emission has the correct spectral shape. 

IRS SH and LH spectral mapping data centered at the posi- 
tion of Vega have been obtained through IRS calibration pro- 
grams (PID 1406, 1409, 1411, and 1413) in 2004. Here we 
present six sets of observations where the slit was placed on 
the Vega position based on PCRS pointing information (no 
IRS peakup). These data were taken in the spectral mapping 
mode with two positions parallel to the slit (7 "47 step" 1 ) and 
three positions perpendicular to the slit direction (l"step _1 ); 
we only used the ones taken at the center perpendicular po- 
sition (two positions alon g the slit center). We used the 
SMART software (v.8.2.5: lHigdon et alJl2004l) to reduce the 
BCD products from the Spitzer Science Center (SSC) pipeline 
version of SI 8. 18. Each of the spectra was extracted with 
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FIG. 1. — IRS spectra of Vega shown in A vs. A format so that a 
Rayleigh-Jeans spectrum is flat. Each observation (AOR) is shifted vertically 
for clear viewing. The unsmoothed, unsealed extracted spectra are shown in 
gray lines, and the scaled, smoothed spectra are shown in colors: red and pur- 
ple for SH while blue and green for LH. The final combined spectrum (black 
dashed line) is also shown in all six AORs for comparison. 

the full-slit mode without sky subtraction, flux calibrated to 
a point source using the SI 8. 18 calibration, and shown in 
Figure Q] We then pinned all the SH spectra to 35.03 Jy at 
10.6 /im accordi ng to the absolute calibration scale defined in 
iRieke et all (120081) and shifted the corresponding LH spectra 
(a scale factor of 0.975) to join smoothly with the SH ones. 
The final spectrum was obtained as a weighted average of 
these six AORs (black dashed line in Figure[T|i. 

The Fomalhaut system was observed with the Spitzer IRS 
instrument in program PID 90 in 2004 June and November. 
Again, we only report here the observations with the IRS 
SH and LH modules, due to saturation of the signal in the 
low-resolution modules. In IRS SH, we obtained a standard 
staring-mode nodded observation, with on-target integration 
time of 503 s. This observation was preceded by a high- 
accuracy pointing peakup on a nearby star with no infrared 
excess, HD 216922, using the IRS Blue (13-19 pm) peakup 
camera. With IRS LH we obtained spectral-mapping observa- 
tions with a strip of nine slit positions separated in the disper- 
sion direction by 4."8 and centered on Fomalhaut. The inte- 
gration time per slit position was 122 s. The spectral mapping 
exercise was preceded by a moderate-accuracy IRS pointing 
peakup using the red (19-26 pm) camera and HD 216922. 
Spectra at the extremes of the spectral map indicated that sky 
emission is negligible compared to that by Fomalhaut, so no 
sky subtraction was performed on the center-position spec- 
trum we discuss here. Data reduction began with basic cal- 
ibrated data products of the IRS data pipeline, version S 1 1 . 
From these data we removed, by interpolation in the spectral 
direction, permanently bad and "rogue" pixels identified in 
the IRS dark-current data for all observing campaigns up to 
and including the one in which each Fomalhaut spectrum was 
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FIG. 2. — Spitzer IRS spectra centered at the star position for Vega (upper panel a) and Fomalhaut (upper panel b). For comparison, Kurucz atmospheric models 
(long dashed lines) and Rayleigh- Jeans spectra (sold red lines) are also shown. The excess spectra (after photospheric subtraction) are shown in the lower panels 
of the plots. The errors for these excess spectra include 2% uncertainty from the photospheric models. 



taken. Again we used SMART for full-slit extractions of spec- 
tra from the two-dimensional data. Similar observations were 
made of a Lac (Al V) in IRS SH and 7 Dra (Kl III). Along 
with template spectra for these stars provided by M. Cohen 
(2004, private communication), we used these observations 
to construct relative spectral response functions (RSRFs) for 
each spectrometer, and in turn used these RSRFs to calibrate 
the Fomalhaut spectra. We estimate the resulting spectropho- 
tometric accuracy to be approximately 5%. 

The final combined, smoothed spectra for both systems are 
shown in Figure [2] For comparison, the stellar photospheric 
models (details see Section I3.ll ) are also shown. Given the 
geometry of the outer cold ring s (radius of 20" in the Fomal- 
haut sys tem (jAcke et aL]l2 012) and radius of 1 1" in the Vega 
system dSibthorpe et al.l l2010)) and the sizes of Point Spread 
Functions (PSFs) in the IRS wavelengths, the spectral flux be- 
yond ^30 /mi is partially contaminated by the cold ring. The 
spectral shape and flux level shortward of ^30 /mi are mostly 
from the star and any unresolved inner warm component. This 
is consistent with the fact that the observed spectrum between 
10 and 13 /im agrees well with the expected stellar photo- 
sphere. The excess spectrum (after stellar photospheric sub- 
traction) is also shown in the lower panel for both systems in 
Figure |2] 

2.2. Herschel PACS Imaging 

Herschel PACS 70 and 160 fim data for Vega and Foma- 
lhaut were obtained by the GT program a nd p ublished in 
an ear ly reduction by Sibthor pe et al.l (12010) and lAcke et al.l 
d2012l) . respectively. We retrieved the archival data and re- 
duced them with the Herschel Inte ractive Pro cessing Envi- 
ronment (HIPE, V9.0 user release, IOt3l2010l) . We applied 
the standard processing steps up to the level 1 stage. Dur- 
ing this process we applied 2nd level deglitching to remove 
outliers with "timeordered" option and 20er threshold. 11 This 
is very effective for data with high levels of coverage. After 
producing level 1 data, we selected the science frames from 
the time-line by applying spacecraft-speed selection criteria 
(between 8"s _1 and 12"s _1 for the slow scan, and 18"s _1 and 
22"s _1 for the median scan). The final level 2 mosaics were 
generated using highpass filtering with the script "photPro- 
ject" and a pixel scale of 1" at 70 /im and 2" at 160 /an. To 

" Details about the timeordered option can be found in the 
HIPS documentation under the PACS photometry data reduction, 
http://herschel.esac. esa.int/hcss-doc-9.0/ 



avoid flux loss in the highpass filtering process, we applied 
a circular mask of 60" radius centered on the position of the 
target. Since the PACS data on Vega and Fomalhaut were 
obtained with different scan rates (10"s _1 and 20"s _1 , respec- 
tively), PSF observations matched to the observing parame- 
ters should be used for comparison. We used PACS data on a 
Boo (Obsld 1342247634 and 1342247635) and a Tau (Obsld 
1342214211 and 1342214212) as our PSF reference and re- 
duced them using the same reduction procedure including the 
masking radius. These two stars are ones of the PACS primary 
calibrators where their fluxes and PSF behaviors are charac- 
terized by Mtill er etaf] (1201 ll) . We have made sure these PSF 
stars have the consistent encircled energy fraction as a func- 
tion of circular ape rture radius derive d by the PACS point- 
source calibration (Mulle ret al.ll2.01 lb . To illustrate the ma- 
jor features seen in the PACS images, we show the final 70 
/im mosaics of Vega and Fomalhaut in Figure [3] along with 
the comparison PSFs. 

3. PHOTOSPHERIC PROPERTIES AND THE IDENTIFICATION OF 
THE CENTRAL UNRESOLVED DISK 

In this section, we first review all ancillary photometry to 
establish the fidelity of photospheric spectrum for Vega and 
Fomalhaut. We then assess whether the excess emission de- 
tected in the IRS spectrum is also detected in the resolved 
Herschel PACS 70 and 160 /im images. Since the stellar pho- 
tospheres dominate the emission in the central part of the re- 
solved images at 70 /im where the systems are best resolved, 
cares must be taken in estimating the contribution of a central 
dust component without involving any further assumption of 
modeling. We do this in two ways: (1) with photometry using 
small apertures that exclude spatially extended emission and 
(2) with PSF subtraction using reference stars to minimize 
the residuals at the star position. The first approach provides 
an estimate of the maximum error on the central component 
due to contamination from other dust emission located in a 
more spatially extended distribution. The second approach 
provides a more accurate estimate of the flux of an unresolved 
source. In both systems, the PACS 160 /im images provide 
only upper limits because the lower angular resolution at this 
longer wavelength makes it difficult to spatially differentiate 
the components. 

3.1. Photospheric Determination Using Ancillary 
Photometry 
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FIG. 3 . — PACS 70 /im images of Vega and Fomalhaut along with their reference PSF stars, a Boo, and a Tau. All images are shown in the same angular scale 
and in the PACS array orientation, i.e., the sub-structures of the PSF are in the same orientation in all four images. The dynamic range of display is from the peak 
value to 1% of the peak value. The color scheme is in logarithmic scale for Vega and a Boo, but in squared root scale for Fomalhaut and a Tau for clarity. In the 
Vega and a Boo images, the two solid circles mark radii of 5" and 7" while the dashed circle marks a radius of 14" (representing the cold planetesimal ring). In 
the Fomalhaut and a Tau images, the solid circles mark a radius of 3". 



The emission of the stellar photosphere from optical to mid- 
IR (^8 /im) was determined in a number of steps. Most 
infrared photometry measurements (like Spitzer/IRAC and 
Akari) are referred to Vega, although Vega is unsatisfactory as 
a standard (i.e., fast-rotating, infrared excess). However, these 
missions also measured Sirius using the same technique, so 
we have taken the measurements with nominal uncertainties 
of 0.01 mag directly compared to Sirius in terms of magnitude 
differences. Sirius (A1V) is very similar in spectral type to 
Vega (A0V) and Fomalhaut (A4V), and is well be haved in the 
infrar ed with no evidence for an infrared excess (tPrice et al.l 
l2004h . 

All three stars are severely saturated in the Two Micron 
All Sky Survey (2MASS) data. Therefore, for accurate mea- 
surements at 2.2 /im, we used data from the DIRBE instru- 
ment on COBE. T he reduction of these data is described by 
iPrice et ail (1201 Oh for analyzing stellar variability. Due to 
the large DIRBE beam (42'x42'), the contribution from stars 
in the field surrounding the target of interest needs to be re- 



moved. We evaluated this effect using 2MASS data, making 
Sirius fainter by 0.014 mag and Fomalhaut fainter by 0.003 
mag while the contribution in the Vega field is negligible. 
Both Vega and Fomal haut are reported to have /T-band ex- 
cess at 1.29%±0. 19% (lAbsil et al.ll2006b and 0.88% ±0.12% 
(lAbsil et al.1 120091) above the photosphere using interferom- 
etry. After accounting for these /T-band excesses of Vega 
and Fomalhaut, the /T-band magnitude difference for photo- 
spheres between Fomalhaut and Sirius is 2.35 mag, and 1.39 
mag between Vega and Sirius. For measurements in the IRAC 
bands, we adopted the results from lMarengo et aLl(l2009l) who 
used PSF fitting technique to recover accurate photometry for 
saturated sources. The magnitude differences between Foma- 
lhaut and Sirius are 2.36 mag, 2.36 mag, 2.36 mag, and 2.33 
mag at 3.6, 4.5, 5.8, and 8 /im, respectively (to be discussed 
further by Espinoza et al. in preparation). The magnitude dif- 
ferences between Vega and Sirius are 1.38 mag, 1.38 mag, 
1.38 mag, and 1.35 mag at 3.6, 4.5, 5.8, and 8 /im, respec- 
tively. Furthermore, Midcourse Space Experiment (MSX ) 
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also measured both Vega and Sirius. The measured 8 /mi flux 
of Vega is lower by 1 % compared to th e MSX predicte d flux 
based on the measurement of Sirius (iPrice et al.M2004b . We 
adopt a magnitude difference of 1.36 mag at 8 /xm between 
Vega and Sirius. 

With the nominal uncertainty for these measurements (0.01 
mag), there is no evidence for excesses from 2.2 to 8 /im (after 
excluding the /T-band excesses from interferometric measure- 
ments) for both systems. We used appropriate Kurucz atmo- 
spheric models for photospheric predictions in other wave- 
lengths by normalizing the model fluxes from 2.2 to 8 /xm to 
the measured photometry that has bee n transferred to th e 
absolute calibration scale proposed by iRieke et all ([2008). 
The parameters in the Kurucz model are T e /f = 9500 K and 
log g =4.0 for Vega, and T eff = 8750 K and log g =4.0 for 
Fomal haut. Given the acc urate Hipparcos parallax measure- 
ments (Ivan Leeuwenl2007l) . the integrated luminosity is 17.45 
Lq for Fomalhaut and 56.05 L Q for Vega (viewing from pole- 
on ). 

Using these normalized Kurucz models, we can then pre- 
dict the photospheric level at the wavelengths of interest and 
determine the excess near the stars by subtracting off the stel- 
lar contribution from the IRS spectra presented in Section l2~T1 
Note that the Kurucz models in the infrared wavelengths are 
basically in the Rayleigh-Jeans regime. 

3.2. PACS Measurements for the Vega Central Source 

The Vega system is viewed pole-on, and its disk has been 
resolved at various wavelengths previously (Holla nd et al. l 
19981; iHeinrichsen et al.1 119981; ISu et all l2005i iMarsh etafl 
2006t ISibthorpe etail l2010t) . As shown in Figure |3] the 
PACS 70 /im image appears to be centrally peaked with a 
smoothed, axis-symmetric extended halo. The size of the 
Vega cold disk observ ed in the higher resolution Herschel data 
( Sibth orpe etai1l2010l) agrees with the one observed in the 
Spitzer data dSu et al.ll2005l) . The stellar photosphere is about 
0.81 Jy at 70 /im, consistent with the centrally peaked mor- 
phology seen at that wavelength. The FWHM of the central 
source is 5 "6x5. "3 measured on a field of 21"x21" centered 
on the star, which is slightly more extended than the mea- 
sured FWHM of the PSF star, a Boo (5."5x5"2). To mini- 
mize the influence of the cold ring (peaked at radius of ~11"- 
14", see Figure|4]l in estimating the photometry of the central 
source, small aperture sizes with appropriate aperture correc- 
tions should be used. On the other hand, the aperture has to be 
large enough to contain most of the flux and to minimize the 
centroid uncertainty. We tried several aperture settings from 
3" to 5 "5 (FWHM) with and without a sky annulus and used 
the PSF star a Boo to derive the values of aperture correc- 
tion. The resultant fluxes range from 1.00 Jy to 1.62 Jy with a 
median value of 1.01 Jy (using an aperture of 5" and sky an- 
nulus between 5"and 7"). At 160 /im, the beam (FWHM) is 
1 l"6x 10." 1 (measured from a Boo), making it very difficult 
to estimate the photometry of the central source alone with- 
out significant contamination from the flux of the cold ring. A 
5" aperture without sky annulus gives a total flux of 0.56 Jy, 
which should be considered as an upper limit since it contains 
some fraction of the cold ring contribution. The photosphere 
of Vega is 0.15 Jy at 160 /im. Therefore, the excess flux at the 
star position, based on small aperture photometry, is 0.2 Jy at 
70 /xm and <0.4 Jy at 160 /im. 

To illustrate our PSF-subtraction approach, cuts were made 
along PA of 120° (corresponding to the horizontal axis of the 
displayed image) of the Vega system with the NW side on 
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FIG. 4. — Surface brightness cuts for the Vega system. These cuts were 
centered at the star position and represent the mean value over a width of 5 
pixels along the PA of 120° (along the horizontal axis in this image). The 
profiles at 70 fj,m are shown in blue color while the ones at 160 (im in red. 
The original cuts (without PSF subtraction) are shown in thick solid lines, and 
the long-dashed lines represent the cuts after the subtraction of photosphere. 
Other thin lines represent the cuts made wi th di fferent levels of subtractions 
(scaling of the PSF). For details see Section |T2"l 

the negative side of the x-axis. These cuts, shown in Figure 
|4] were centered at the stellar position at 70 /im, and repre- 
sented the mean value over a width of 5 pixels (i.e., 5" and 
10" at 70 and 160 /im, respectively). No significant asym- 
metry in disk brightness nor center offset were seen in these 
cuts. After subtracting a photospheric scaled a Boo PSF, an 
additional point-like source is clearly present in the resultant 
cut (blue dashed line in Figure 2]i at 70 /im. The maximum 
of the point-source contribution (star and unresolved disk) can 
be estimated by forcing the peak values (single pixel) matched 
between the Vega and a Boo data. A flux density of 1.19 Jy 
for such a PSF subtraction is required and its resultant profile 
is shown as the thin solid blue line in Figure [4] suggesting a 
maximum brightness of <0.38 Jy for this unresolved disk at 
70 /im. Using the flux obtained with the small aperture pho- 
tometry (0.2 Jy for the unresolved disk), the resultant profile 
is shown as blue dotted-dashed line in FigureH] relatively flat 
in the central ±5" region as expected. Therefore, the bright- 
ness of the unresolved disk in the Vega system is ^0.2 Jy 
with a maximum value <0.38 at 70 /im. Du e to the uncer- 
tainty of the PSF (as much as 10% at 70 /xm; iKennedy et alJ 
1201 2l) and that of flux calibration, we simply assume a conser- 
vative lower-bound error of 20% (i.e., a minimum flux of 0. 16 
Jy at 70 /im for the central unresolved disk). At 160 /im, the 
maximum scaled PSF is ^0.45 Jy by forcing the peak value 
to zero, suggesting a maximum brightness of <0.3 Jy for this 
unresolved disk component. 

3.3. PACS Measurements for the Fomalhaut Central Source 

Since the Fomalhaut debris system is inclined by 67°, it is 
very difficult to separate the point source from the cold ring 
even using very small apertures. One can still estimate the 
maximum brightness of the unresolved component by assum- 
ing there is no contamination from other sources inside a very 
small aperture. For an aperture radius of 3 ", the encircled flux 
is 0.95 Jy at 70 /xm and 0.77 Jy at 160 /tm after applying ap- 
propriate aperture corrections. Taking out the contribution of 
the stellar photosphere (0.368 Jy and 0.07 Jy at 70 and 160 
/xm, respectively, based on our stellar photospheric model), 
the unresolved disk component has a maximum flux of <0.58 
Jy and <0.7 Jy at 70 and 160 /xm, respectively. 
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FIG. 5. — Surface brightness cuts along the major axis of the Fomalhaut 
system. These cuts were centered at the star position (marked as the solid 
vertical line) and represent the mean value over a width of 4 pixels along the 
minor-axis. Line styles and colors are similar to those used in Figure [4] The 
dotted and dotted-dash vertical lines mark the peak positions of the ansae at 
70 and 160 /im, respectively, while the long-dashed vertical line represents 
the mid-point of the cold ring. 

Similar to Vega, cuts were made along the major-axis of 
the disk and shown in Figure [5] These cuts were centered at 
the stellar position at 70 /im, and represented the mean value 
over a width of 4 pixels (i.e., 4" and 8" at 70 and 160 /im, 
respectively). The center of the outer cold ring (the dashed, 
vertical line in Figure [5]l was estimated by the mid-point of 
the two bright peaks (marked as dotted, vertical lines for 70 
/im and dot-dashed, vertical lines for 160 /im in Figure [5]l at 
both wavelengths. The 160 /im ansae peak closer to the ring 
center than the o nes at 70 /im, which was first reported by 
lAcke et ail d2012f) and suggested to be due to blurring in the 
large PSF at 160 /im. Although the peak positions at 70 and 
160 /im are not the same, the cold ring centers at the same po- 
sition at both 70 and 160 /im. We used a PSF star, a Tau, for 
photospheric subtraction by scaling it to match the expected 
photospheres (0.368 Jy and 0.07 Jy at 70 and 160 /im). The 
photosphere-subtracted cuts are shown in long-dashed lines 
in Figure [5] At 70 /im, it is clear that an additional source 
of emission is present at the center of the disk. We tried two 
different PSF subtractions to estimate the brightness of this 
additional component. First, we arbitrarily increased the scal- 
ing of a Tau, fixed in the stellar photosphere position, until 
the central region (from -5" to +5" ) has a relatively fiat dis- 
tribution in the cut. The resultant cut is shown as a thin solid 
blue line with an additional scale of 0.17 Jy. Second, we sub- 
tract a second a Tau PSF by adjusting both the position and 
the scaling until the resultant cut is relatively fiat (dot-dashed 
line in Figure [5J. The second method gives a scale of 0.165 
Jy for this additional source. Combining both methods, we 
conclude that the central unresolved component is 0.17 Jy at 
70 /im. There is no easy way to estimate the error associated 
with this number since it depends on the detailed structures of 
various components, we simply assume a 20% error for fur- 
ther analysis. At 160 /im, it is not obvious that an additional 
source is required in the system. We estimated the upper limit 
of such a component by increasing the scaling factor of a Tau 
fixed at the star position (any potential offset between the star 
and this inner component is washed out by the large beam size 
at 160 /im). An upper limit of 0.08 Jy at 160 /im was inferred. 

In summary, the central unresolved disk component in the 
Fomalhaut system is about 0.17 Jy at 70 /im with upper- and 
lower-bound fluxes of 0.58 Jy and 0.136 Jy (assuming 20% 



uncertainty), and <0.08 Jy at 160 /im, based on the analysis of 
the data only (without invoking any assumption of modeling). 
Our values are c onsistent with the best-fit model presented in 
lAcke et all d2012l) . Their model estimates that the fluxes of 
the central point source (star + unresolved disk component) 
are 0.54 Jy (5% of the total flux) and 0.124 Jy (2% of the 
total flux) at 70 and 160 /im, respectively, implying that the 
unresolved disk accounts for flux of 0.172 Jy at 70 /im and 
0.054 Jy at 160 /im. 

4. ANALYSIS 
4.1. Unresolved Warm Excess in the Vega System 

The Vega Spitzer 24 /im observation dSu et al.ll2005l) was 
severely saturated, making it difficult to constrain the bright- 
ness of this unresolved component without invoking some 
modeling assumptions. Therefore, we seek other relevant 
measurements in the mid- infrared to valida te the excess levels 
seen in the IRS spectrum. Tokunagal d 19841) measured a hand- 
ful of nearby A-type stars at 20 /im using the IRTF bolometer 
with a beam size of 5" to define the 20 /im magnitude sys- 
tem (relative to Vega). Only two stars, a CMa (Sirius) and 
7 UMa, in his list do not have a 24 /im excess. We used the 
color V-[20] of these two stars to extract the excess of Vega 
at 20 /im. We adopt V of -1.40 and 2.40 mag for a CMa and 
7 UMa, suggesting a photosphere color V-[20] of -0.04 mag 
in Tokunaga's system. Based on the observed V-[20] color of 
0.03, the excess of Vega is ^7% above the photosphere (i.e., 
0.7 Jy). It is difficult to assess the errors associated with the 
20 /im measurement; therefore, we simply assume a maxi- 
mum 50% error at this wavelength. Another source of mid-IR 
measurements for Vega comes from MSX photometry where 
the excesses (rel ative to Sirius) at 14.65 and 21.34 /im have 
been reported by iPrice etail (120041) to be 4% and 17% above 
the photosphere (i.e., excesses of 0.7 Jy at 14.65 /im and 1.5 
Jy at 21.34 /im), respectively. The contamination from the 
cold ring needs to be taken into account given the large beam 
size of MSX (a resolution of ^20")- Assuming a typical dust 
temperature of 60 K for the cold ring and normalizing the cold 
ring flux observed in the far-infrared dSu et al.ll2.005b , the flux 
contamination from the cold ring is less than 0. 1 % of the pho- 
tosphere at 14.65 /im, and ^4% of the photosphere at 21.34 
/zm. The nominal error for the MSX measurements is ~1,5% 
( Riek e et al.l2008l) . which includes the errors from stellar pho- 
tospheric prediction. Including the possible contamination 
from the cold ring, the final errors are 0.31 Jy and 0.38 Jy at 
1-cr for the excesses at 14.65 and 21.34 /im. The Vega system 
was rece ntly observed by the Submillimeter Array (SMA) at 
880 /im (Hugh es et al.ll2012l) . Interestingly, the spatially un- 
resolved 880 /im flux within 5" of the Vega photosphere lies 
slightly above the predicted photospheric value (see Figure|6]), 
although only at the less than 2a level. Combining with our 
estimates from the PACS images and the IRS spectrum, the 
excess SED of the central compo nent in the Veg a system is 
shown in Figure [6] We note that lLiu et al.l d?004) also report 
an upper limit (^250 mJy) of 0.7% (la) of the photospheric 
level at 10.6 /im, using nulling interferometry (BLINC/MMT) 
that probes a region less than 1 "5 (12 AU) from the star. 

The spectral shape (A <30 /im) of the inner component is 
well presented by a blackbody of 170 K, just above the tem- 
perature at which water ice sublimates in vacuum. This tem- 
perature corresponds to a distance (^2.7 AU) in the middle of 
our asteroid belt in the solar system, but ^14 AU in the Vega 
system (using an average stellar luminosity of 37 L Pl viewed 
by the dust along the equator of the star; lAufdenberg et al.l 
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FIG. 6. — Excess SED of the inner warm component in Vega. The photo- 
spheric subtracted IRS spectrum is shown as green dots with errors (including 
the 2% errors from the photosphere) shown as gray area; other IR and sub- 
millimeter excess photometry are shown as blue squares. For comparison, 
blackbody curves of 170 K are also shown with B\ as the long-dashed line 
and A-°- 5 B A as the dotted-dash line. The 10.6 /Urn BLINC/MMT \cr upper 
limit is also shown. 

2006) assuming blackbody-like emitters. 

4.2. Unresolved Warm Excess in the Fomalhaut System 

The existence of a close (unresolved) warm compo- 
nent around the Fomalhaut disk was first suggested by the 
Spitzer 24 observation where an additional point-like 
source (0.6±0.2 Jy) along with the expected stellar pho- 
tosp here is required to fit the resolved disk image at 24 
/im (jStapelfeldt et al. 2004). The central component of the 
Fomalhaut system was als o tentatively de tected by the ALMA 
observation at 850 /im (iBoley et al.ll2012h . A flux density of 
~4.4 mJy at the star position was estimated after the primary 
beam correction, which is quite uncertain at the star position 
because the ALMA observation was centered at Fomalhaut b. 
The expect ed photosphere at 850 fim is ^2 mJy (not 3 mJy 
as stated in iBolev et alj|2012). Therefore, the excess at 850 
/im could be as high as 2.4 mJy. We simply took this at a 
face value, and assumed a 50% error (as maximal). The ex- 
cess SED for the Fomalhaut central unresolved component is 
shown in Figure [7] along with the excess spectrum measured 
by IRS. 

We note that lAcke et alJ (1201 2l) suggested ionized gas (free- 
free emission from a hot stellar wind) being responsible for 
the unresolved central excess from K band to 70 /im, and 
derived F v oc ^ 8±0 1 power law. With their derived power- 
law, the excess at 850 /im (350 GHz) is expected to be 18-23 
mJy, which is clearly not consistent with the ALMA observa- 
tion. In addition, this power-law spectrum is inco nsistent with 
the 2.2-5.8 /im colors discussed in Section [37X1 Instead, we 
suggest that the excess emission detected longward of ~13 
/im arises from thermal dust emission. Similar to the warm 
component in the Vega system, the spectral shape (A <30 /im) 
of the inner component in Fomalhaut is also well representa- 
tive by a blackbody of 170 K, suggesting a radial distance of 
~1 1 AU from the star for blackbody-like grains. 

4.3. Asteroid-belt Analogs 

The warm excesses detected in both systems are restricted 
to the vicinity of the star; i.e., unresolved at various wave- 
lengths. The amount of warm excess emission is derived ei- 
ther from resolved images where the central component is 
(mostly) separated from the cold ring or from the spectra 
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FIG. 7. — Excess SED of the inner warm component in Fomalhaut. Sym- 
bols and lines are similar to what were used in Figure [6] The upper error 
bar at 70 fim represents the (unlikely) maximum value estimated from small 
aperture photometry and the lower error bar is 20% lower than the estimate 
from PSF subtraction. Other error bars are shown as l-c For comparison, 
blackbody curves of 170 K are also shown with B\ as the long-dashed line 
and \~ 3 B\ as the dotted-dash line. 

taken centered at the star position. Given the sizes of the IRS 
slits, some amount of the flux from the cold rings is expected 
to fall into the slit, especially at wavelengths longer than ^30 
/im. It is difficult to estimate the exact flux contamination 
without further modeling because it depends on the geome- 
try of the cold ring at different wavelengths. One can esti- 
mate the maximum flux contamination in the worst scenario 
case by assuming the warm and cold components are spatially 
coincident. Therefore, the maximum flux contamination is 
<1% of the photosphere for A <20 /im and could be up to 
50% at 33 /im given the typical dust temperature of the cold 
ring (^70-50 K). Thus the derived dust temperature (~170 K 
in both systems) based on the data shortward of 30 /im has 
very little contamination from the cold ring. In other words, 
the excess arises from material less than 3" (unresolved with 
an FWHM~6") from the star in both cases. Similar to the 
definition of habitable zones around stars, debris disk struc- 
tures should be identified and characterized in terms of dust 
temperatures rather than physical distances so that the heating 
power of different spectral type of stars is taken into account 
and common features in disks can be discussed and compared 
directly. The characteristic temperature of ^170 K for the ex- 
cesses suggests they are asteroid-belt analogs. It corresponds 
to a distance of 2"(~14 AU) from Vega (using the lower lu- 
minosity viewed by the dust in the equatorial direction) or of 
li'S (~11 AU) from Fomalhaut assuming blackbody radia- 
tors. The resultant location of the asteroid belt given a dust 
temperature depend s greatly on grain pro perties. Using as- 
tronomical silicates (Laor & Drainelll993l) . the dust can be as 
close as l."3 (10 AU) from Vega, and 1" (8 AU) from Fo- 
malhaut using grains with a radius of 10 /im. Furthermore , 
the IRAM PdBI observations at 1.3 mm bv lPietu et aTJ d201 ll) 
do not find any exc ess within 2"5 while the SMA data by 
iHuphes et alJ d2012l) suggest a tentative excess within 5" from 
the star, suggesting that the warm excess around Vega likely 
arises from emission outside a radius of l."3 (10 AU) from 
the star. The fractional luminosity is 7xl0~ 6 for the Vega 
warm component and 2 x 10~ 5 for the Fomalhaut warm com- 
ponent, muc h more luminous tha n our current zodiacal cloud 
(TO" 8 - lQ- 7 : lDermott et ai]l2002h . 

The origin of the warm excess in the vicinity of the star is 
unknown. One hypothesis is that the warm excess originates 
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from dust grains in the cold belt, which are transported in- 
ward by P-R and/or stellar wind drag, a s suggested for the 
outer w arm belt of e Eri debris system by Reidemeiste r et all 
d201 lb . In the Fomalhaut system, the presence of an in- 
ner disk, extending inward up to —35 AU (4"5) and com- 
posed of grains dragg ed-in from the co ld planetesimal belt, 
has been suggested bv lAcke et al.l (12012) to explain the PACS 
70 /im profile. Conventionally, a PR- transported flow has a 
constant surface density and is expected to extend to the star, 
resulting in a surface brightness profile that is centrally peaked 
at the star. However. lAcke et aD,(l2012l) find that the inner disk 
in Fomalhaut has a clear truncation (outside the unresolved 
warm component). Therefore, neither the warm nor the hot 
excess observed in the Fomalhaut system is unlikely to have 
arisen from the inward transport of grains by P-R drag. On the 
other hand, we cannot rule out such a possibility (either solely 
or partially due to the dragged-in grains) in the Vega system 
(it is centrally peaked at 70 /im) without detailed modeling 
for the whole system at multiple wavelengths (K. Y. L. Su et 
al. in preparation). Nevertheless, the observed dust tempera- 
ture (170 K) of the warm excess suggests this component does 
not extend all the way to the star; instead, an inner truncation 
(>1") is required. 

Furthermore, the tentative detection of the central com- 
ponent at submillimeter wavelengths in both systems argues 
for the presence of large grains, suggesting the warm emis- 
sion likely arises from a full spectrum of particle sizes sim- 
ilar to the cold planetesimal belt, i.e., a form of planetes- 
imal belt like our asteroid belt. The wavelength-dependent 
power index at long wavelengths (F v ~ A"') is a measure of 
the grain size distribution in a collision-d ominated debris disk 
dWvatt & Denj|2002tlGlspar et al.ll2012l) . The cold planetes- 
imal belts observed in bright debris disks have much steeper 
slope s compared to Rayleigh- Jeans slope (e.g., iGaspar et al.1 
2012). Similar behavior is seen in the warm component in the 
Vega and Fomalhaut systems (see Figures [6] and |7J, implying 
the presence of large grains. Limited by the uncertainties in 
the current observations, the exact slope of the warm com- 
ponent at long wavelengths cannot be determined accurately. 
Future high-resolution observations in the submillimeter and 
millimeter will help shed light on the nature of the warm ex- 
cesses. 

One interesting note on the nature of these warm excesses 
is that the observed levels of dust (fd ~7x 10~ 6 for Vega and 
/j~2x 10~ 5 for Fomalhaut) are consistent with being the 
in situ, steady-state collisional evolution of large parent bod- 
ies for the lifetime of the stars (—400 Myr). The maximum 
fraction luminosity (/ mflA ) for such a sy s tem c an be estimated 
using Equation (18) from IWvatt et al.l d2007l) where f max is 
on the order of (2-3) x 10~ 5 for both systems assuming an as- 
teroid belt at —10 AU with a width of — 1 AU and consisted 
of — 10 km (diameter) planetesimals in collisional cascades. 
Due to the uncert ainties of some parameters in the model, 
IWyatt et aD (l2007h suggest that a transient event producing 
the observed dust is only required when f„b s 3> I0 3 f max . Thus, 
the observed dust in the warm component is consistent with 
it being generated through collisional grinding in an asteroid 
belt in both systems. 

5. DISCUSSION 

Vega and Fomalhaut really live up to their names as debris 
disk twins: both possess a hot 2 /im excess revealed through 
interferometric observations, a cold (—50 K) belt analogous 
to our Kuiper belt that has been in the spotlight of space in- 



frared facilities, and a warm (—170 K) belt analogous to our 
asteroid belt. This warm-and-cold- belt debris arch itecture has 
also been identified i n HR 8799 (Su et al. 2009) and e Eri 
dBackman et al.l 120091) . In Section IBTTl we first discuss the 
possible origin of these two-belt systems in terms of their for- 
mation and evolution. Taken together with the known proper- 
ties of planets and mass limits from non-detection, we spec- 
ulate that the large gap between the two belts is likely to be 
sustained by multiple (low-mass) pl anets . We then review the 
properties of 2 /im excess in Section lBTZl and suggest an alter- 
native explanation for the hot dust component. 

5.1. Origin of the Two-belt Systems 

In our solar system, the locations of the minor bodies that 
failed to form planets are elegantly arranged and sculpted 
by planetary perturbations over the course of 4.5 Gyr evo- 
lution. The inner edge of the Kuiper belt's dusty disk is 
thoug ht to be maintained by massive planets dLiou & Zookl 
1999), whereas the more tenuous asteroid-belt dust (i.e., zo- 
diacal cloud) has a structure also determined by gravita- 
tional perturbations via both the giant an d terrestrial planets 
dDermottet al.lll994t iMurray et al1ll998l) . It has been sug- 
gested that the dominant source of dusty debris inside —5 
AU in ou r solar system result s from the breakup of aster- 
oids (e.g., Dermott et al. 2002). However, a recent dynami- 
cal model incorporating multip le sources (ast e roid s and short- 
and long-period comets) by iNesvornv et al.l (1201 Ol) suggests 
that particles originated from Jupiter-family (short-period) 
comets dominate the mid-infrared emission in the zodiacal 
cloud while the contribution of asteroidal dust is <10%. The 
origin of our own zodiacal cloud is still a matter of consider- 
able debate, making understanding of the warm dust around 
nearby stars (i.e., exo-zodi) even more valuable. 

Dust locations solely based on temperatures derived from 
ex-solar debris disk SEDs are ambiguous. However, based 
on resolved images at multiple wavelengths, four systems - 
Vega, Fomalhaut, e Eri, 12 and HR 8799 - clearly have two 
dust belts. The disk SEDs of these four systems are shown 
togethe r in Figure [8] for easy comparison. As discussed in 
Section 14.31 the dust observed in the warm components for 
both Vega and Fomalhaut is consistent with being generated 
in a steady-state collisional cascade in a planetesimal belt. A 
similar conclusion obtains for the other two systems. For HR 
8799, the warm dust component has a t emperature — 150 K 
and fractional luminosity f d — 2 x 10" 5 dSu et al.ll2009h : the 
latter is much less than the maximum fractional luminosity 
for steady state collisional dust production, f max w3x 10 , 
estimated for a belt at 8 AU at HR 8799's age of -30 Myr. 
For the e Eri system, the war m component has fd — 3 x 10~ 5 
and temperature of —150 K (Backman et al. 2009); the latter 
is compatible with f max ps 8 x 10~ 6 for a belt at 3 AU at an age 
of —800 Myr. Based on the similarity of the observed SEDs 
and resolved disk structures, it is very likely that the warm 
components in these four systems originate in well-separated 
planetesimal belts (remnants of planet formation), resembling 
the configuration of our own solar system which has two left- 
over planetesimal belts: an asteroid belt near the water-frost 
line, between Mars and Jupiter, with a characteristic temper- 
ature of -170-150 K, and a Kuiper belt at -30-55 AU with 
dust emission peaked at —50 K. 

12 For the e Eri disk, we refer the inner warm belt that dominates the 
emi ssion at 24 (im rather the outer warm belt that emits prominently at 55 
(Backman et al. 2009). 
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FIG. 8. — Excess SEDs for e Eri, Fomalhaut, Vega and HR 8799 where 
measurements are shown in various symbols with colors differentiating the 
objects and large size symbols being the integrated photometry while smaller 
ones representing spectroscopic measurements. The photometry of the warm 
components is also shown as the smaller symbol. The excess of each system 
is well represented by two (warm as long-dashed line and cold as dashed line) 
blackbody emissions: e Eri: 150 K and 50 K, Fomalhaut: 170 K and 50 K, 
Vega: 170 K and 50 K, and HR 8799: 150 K and 45 K. Note that the cold 
component is modeled with modified blackbody as \~ 9 B\ in order to fit the 
sub-millimeter observations. 

There are two aspects of the two-belt systems that must 
be explained. The first one is related to how these two sep- 
arate belts were created, and the other one addresses how 
the system maintains such a large gap without dust filling 
in it. Mechanisms to explain both the formation and evolu- 
tion of the two belt systems are required to fully account for 
the observed pattern. Planetesimal belts can be expected in 
regions where there was not enough material or not enough 
time to form a planet, because once formed, a planet would 
scatter or accrete the surrounding planetesimals. Therefore, 
the stable locations of left-over planetesimals are governed 
by where the giant planets form and their migration history, 
which may include strong dynamical instabilities. Conse- 
quently, the location of the observed dust is not expected to 
show much order among systems since it greatly depends on 
the numbers and positions of the giant planets and their mi- 
gration history. In addition, it is an observational fact that 
higher mass stars harbor a more massive pr otoplanetary disk 
dNatta et all 120001; IW illiams & Ciezal l201ll) , and more giant 
planets dJohnson et al.ll2.01Qh . Generally, one does not expect 
to find similarity in planetary configurations between low- 
mass (i.e., low luminosity) and high-mass (i.e., high luminos- 
ity) stars. 



Among the four systems shown in Figure [8] only e Eri is 
a solar-like star and the rest are early-type stars. The lumi- 
nosity (the dust heating power) difference ranges roughly two 
orders of magnitude among them, and yet the characteristic 
dust temperatures of the warm and cold belts are very simi- 
lar. It is interesting to note that many unresolved systems also 
have a similar two-belt configuration in the dust distribution 
based on t heir temperatures derived from unre solved excess 
emission (IChen et al.l2009tlMorales et al.l201 ll) . The fact that 
the temperature of the warm belts peaks at similar temper- 
atures (~ 170-1 90 K) bet ween early-type and solar-like sys- 
tems (Mor ales et alJ l20TTh . however, suggests temperature- 
sensitive mechanisms play a major role in determining plan- 
etesimal configuration in the warm belts. Giant planets are 
expected to form right outside the water-frost line where the 
amount of solid mat erial and dyn amical timescales favor the 
formation process (iKretke & Linll2007h . 

The young HR 8799 system has four giant planets 
(iMarois et al. 2008, 2010) separating the inner and outer dust 
belts dSu et alj r2009). These four planets are very massive 
(<~7— 10 Mj) so that they dominate the dynamics in the HR 
8799 system. Dust particles spiraling inwards from the cold 
belt under P-R drag are likely to be dynamically scattered and 
ejected by one of the planets before they reach the inner re- 
gion. Thus, the warm component in the HR 8799 system 
is unlikely to arise from grains generated in the cold plan- 
etesimal belt. Comets, on plunging orbits originating from 
the cold belt could possibly cross over and then break-up or 
sublimate to deliver dust interior to the giant planets. How- 
ever, it is difficult to estimate the inward flux of cometary 
bodies inside the orbits of the giant planets with out detailed 
numer ical simulations such as the one done by iBonsor et al.l 
(|2012|). Furthermore, the comet delivery scenario proposed 
by lNesvorny et al.l ([20 1 ) might work for our solar system at 
its current age, it is not clear that this scenario is directly ap- 
plicable for younger systems (such as the ones we discussed 
here) where in-situ dust generation in the younger, more mas- 
sive planetesimal belts is expected to dominate. 

In the remaining three systems, Fomalhaut is the only 
ot her one that harbo rs a directly detected planet, Fomalhaut 
b ( Kal as et al] l2008). Although the reality of Fomalhaut b has 
been questioned because the spectrophotometry of Fomalhaut 
b cannot be reconciled with models for thermal emission from 
a young giant pl anet, a problem originally acknowledged in 
iKalas et all d2008l) and further d e monstrated by Spit zer non- 
detections (Mareng o et al.ll2009t iJanson et al.ll2012l) . a per- 
turbing plan et is required to a c count for the o bserved disk 
asymmetry (IKalas et alJ 12001 lOuillenl l2006t IChiang et alJ 
I2009I) . Nevertheless, rece nt r e-analyses on the p ublic HST 
databy lCurrie et al.l(l2012l) and lGalicher et al.l (120121) confirm 
the detection of Fomalhaut b with its position found to be just 
interior to the ring and comoving with the star, making it a 
candidate to maintain the sharp inner boundary of the Foma- 
lhaut cold ring. The orbit of Fomalhaut b has not yet been 
fully demonstrated to be consistent with the range of orbits 
required for the perturber, and its mass is also quite uncer- 
tain due to lack of detections at other wavelengths, ranging 
from <a few Mj (from the constraint of cold belt as detected 
in scattered light, Chiang et al. 20091) to a few M ffi (from the 
constraint deri ved from the pro perties of the cold belt in the 
submillimeter. lBoley et al.l2012l) . Furthermore, ground-based 
high-contrast observations have also s et a mass limit of <2 M 



for any planet between ^10^-0 AU (Kenworthv et al. 
and of < 10-16 Mj from 3 to 10 AU (IKenworfhy et al 



2009) 



2013) 
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in the Fomalhaut system. 

The clumpy structures claimed from millimeter and sub- 
millimeter imaging of the Vega cold Kuiper-belt analog were 
often ta ken as signatures of gravitational per t urbations by 
nlanets (iHolfand et alJ[T998t IWilneret alJl200l IWv^lIOOl 
Marsh ^taTI 120061) . Recent observations at higher resolu- 
tion and sensitivity fail to detect the clumps, and instead 
are consistent wi t h a smooth, broad, and axisymmetric disk 
dPietu et al.H20Tll:lHughes et al.ll2012h . In any case, Vega has 
been the target of several deep searches for planets through di- 
rect imaging, and strong limits have been placed at H, L', and 
M bands f or planets with masses > 3 -4 My betwe en —20 and 
-70 AU dMarois et al.ll2006t iHeinze et alJl2008l) . For e Eri, 
planets with masses >2-3 My in radia l distance of 6-35 AU 
are ruled out at H , L', and M bands (lLafreniere et al.l 120071: 
Heinze et all 120081) and Spitzer IRAC bands dMarengo et all 
2006, 2009). Using radial velocity (RV) and astronometry 
techniques, a close, eccentric planet, e Eri b, with a semi- 
major dista nce of —3-4 AU wa s identified by Hatz es et"aT] 
(2000) and iBenedict et all (120061) : although the discovery of 
the inner warm dust belt (—3 AU) casts doubt on th e existence 
of a high eccentricity planet dBackman et al.ll2009l) . Further- 
more, a recent analysis of all av ailable RV data of e Eri by 
Anglada-E scude & Butlerl (1201 2|) found a significantly differ- 
ent orbital solution and suggested that the long-term RV vari- 
ability is likely due to stellar activity cycles rather than a pu- 
tative planet. Nonetheless, e Eri b, if it exists just outside 
the warm belt, is a ca ndidate to shepherd the warm dust belt 
dBackman et alJl2009h . 

Although the exact locations of the warm components 
in these disks are unknown (unresolved), the orbital ratios 
(Rcoid I Rwarm) between the outer cold belt (mostly resolved in 
the submillimeter and millimeter wavelengths) and the warm 
belt estimated from the dust temperature are roughly >10 in 
all four systems, 13 similar to the ratio in our own solar sys- 
tem (the asteroid belt at —3 AU and the Kuiper belt at —35 
AU). A significant deficit of large grains (best tracers for the 
unseen parent bodies) in the region between the belts are evi- 
dent due to the fact that high-resolution submilli meter obser- 
vation s do not detect such a filled-in component (Boley et al. 
12012 ), although sma ll grains can drift inward from the cold 
belt dAcke et"al]|2012l) . The zone between the warm and cold 
belts that is mostly free of dust is very large. If the observed 
dust is being generated in both belts through collisional cas- 
cades of large parent bodies, we need a cleaning mechanism 
to maintain such a large dust-free zone. 

The large gap between the warm and cold belts may be 
maintained by one or multiple planet-mass objects in the gap. 
If these planets have dynamical influence over the entire gap, 
they are likely to scatter any material that is either generated in 
the gap or drifts into it. The dynamical influence of a planet is 
given approximately by the extent of the overlap of first-order 
resonances which creates an un stable "chaotic zone" in the 
vicini ty of a planetary pertur ber (lWisdomll980l iDuncan et al.l 
1989t lMustill & Wyatfl20T2h . Applying this criterion, we find 
that in the Vega system, the large gap could be maintained by 
a single hypothetical object with mass of a few 100 My (i.e., a 
star/brown dwarf) in a circular orbit; for Fomalhaut, a similar 
mass estimate of a few 100 My obtains for a single hypotheti- 
cal object in an eccentric orbit (e = 0.1, based on the observed 

13 Vega: R warm ~14 AU and R coU -110 AU; Fomalhaut: R„ mm -10 AU 
and R coU ~140 AU; HR 8799: R war „, ~10 AU and R coM ~100 AU; t Eri: 
Rwarm ~3 AU and R coU ~35 AU. 
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FIG. 9. — Mass-dependent chaotic zones for two equal-mass planets in 
the Vega (left panel (a)) and Fomalhaut (right panel (b)) systems, related to 
the boundaries of the inner and outer belts (shown as black, dashed, vertical 
lines). In both panels, the orbital radius of the inner planet is shown as red 
solid lines while the one of the outer planet is shown as blue solid lines. In 
Vega, the chaotic zone width is computed using e =0 for the planet's eccen- 
tricity, but for Fomalhaut we adopt e =0. 1 . The dotted-dash lines represent 
the outer boundary of the chaotic zone for the inner planet while the long- 
dashed lines represent the inner boundary of the chaotic zone for the outer 
planet. For two equal-mass planets, the planet-star mass ratio, fj,, has to be 
greater than 0.015 (0.03) for these two planets to maintain the dust-free zone 
in the Vega (Fomalhaut) system. The current observational upper limit for 
the planet mass is shown as the horizontal dotted lines. 

eccentricity of its cold dust ring). One can also imagine a 
lower mass perturber in very eccentric orbit being responsi- 
ble for maintaining such a large gap; the required eccentricity 
can be estimated by assuming that its pericenter and apocenter 
are near the locations of the warm and cold belts. We find that 
a large eccentricity (e —0.8) is required, implying an object 
with a mass of —50 My in both cases. Our simple estimates 
show that the gaps in both Vega and Fomalhaut are too large 
to be explained by a single circular/eccentric perturber, with- 
out contradicting the current upper limits of a few My based 
on non-detections of planetary companions in the gaps. 

The minimum number of planets residing between the two 
belts is two - one inner planet outside the warm dust belt to 
shepherd the inner planetesimal belt, and one outer planet in- 
terior to the cold dust belt to scatter large grains that drift 
inwards from the cold belt. We can then estimate the mass 
of two equal-mass planets in this case. This is illustrated in 
Figure [9] where we plot the orbital distances of the two equal- 
mass planets and their associated chaotic zone widths as a 
function of the mass ratio between the planet and the host 
star. In the case of Vega (M» = 2.5M Q and we assume e = 
for the hypothetical planets), we find that two —40 My plan- 
ets are needed to maintain the dust free zone. In the case 
of Fomalhaut (M„ = 2.0M Q and we assume e = 0.1 for the 
hypothetical planets), we find that two —63 My planets are 
needed. Although these masses are likely overestimates be- 
cause the single-planet chaotic zone formulae do not account 
for the strong secular perturbations that ca n extend the unsta- 
ble zo nes in multiple planet systems (e.g.. lMoro-Martm et al.l 
2010), it is evident that just two planets in low eccentricity or- 
bits having mass — My are inadequate for explaining the gaps. 

Our simplified approach, in estimating the masses for sin- 
gle planetary perturbers and for two equal-mass perturbers for 
explaining the large gaps in Vega and Fomalhaut, gives lim- 
its that are more than an order of magnitude higher than the 
—few My planet mass limits for these systems based on non- 
detections of planetary companions in the gaps. It is, there- 
fore, probable that multiple, lower mass planets are responsi- 
ble for clearing the gaps. We can estimate the numbers and 
orbital radii of multiple, Jupiter-mass planets that would clear 
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this region, by assuming that they are separated such that their 
chaotic zones just fill the regions between them. We find that 
four or five IMj mass planets are required for the large gap in 
Vega and in Fomalhaut. (This is reminiscent of the HR 8799 
system in which four planets separate the two dust belts.) Up- 
coming ground-based high-contrast, direct imaging surveys 
using facilities like LBT/LBTI, Gemini/GPI, VLT/SPHERE 
should be able to find and/or place tighter mass limits on the 
planets in the large gap. Stars that have a debris disk with 
two distinct dust belts separated by a large gap are attractive 
targets for future searches of exoplanets. 

5.2. The Mystery of the Hot Dust 

The presence of hot dust in the close vicinity of both sys- 
tems is revealed by ground-based interferometric observa- 
tions. For Fomalhaut, 0.88% ±0.12% of excess emission in 
the K band was reported by lAbsil et al.l (120091). For Vega . 



2006) 



2011) 



a 1.29%±0.19% excess in the K band (lAbsil et al.l 
and 1.23% ±0.45% excess in the H band dDefrere et al 
were reported. Generally, these studies have ruled out other 
possible sources such as low-mass comp anions or stellar 
winds to be the cause of excess emission. Men nesson et alj 
(1201 ll) also put a very tight constraint on the location of 
the hot dust (withi n 0.2 AU) by combin ing all interferomet- 
ric measurements. Defr ere et al.l (1201 lb further modeled the 
properties of the hot excess of Vega (both spatial visibility 
and SED) and reached several conclusions: the dust responsi- 
ble for this hot excess (1) has very steep density and particle 
size distributions (density power index less than -3 and par- 
ticle size power index ~ -5), (2) resides from 0.1 AU (dust 
sublimation radius with temperature of 1700 K) to less than 
0.2 AU from the star, and (3) has a minimum particle size of 
0.01-0.2 /mi with a significant fraction of carbonaceous com- 
position (due to the fact that carbonaceous grains have a much 
higher sublimation temperature and lack prominent miner- 
alogical features). Although their model does include large 
grains (with maximum size of 1000 /im) in the calculation, 
there is basically no grain with sizes larger than 0.3 /im (the 
average size particle is 0.27 /im) with such a steep size dis- 
tribution. In other words, the only explanation for such a hot 
excess is a population of sub-/im carbonaceous grains located 
in a narrow -ring-like region at the dust sublimation radius (0. 1 
AU). In the case of Vega, this hot excess is very different from 
what we know about the zodiacal cloud in our solar system, 
which has a rel atively flat density d istribution (density power 
index of -0.34; lKelsall et alJI 19981) a nd mostly contains larg e 
particles with sizes of ^10-100 /im dFixsen & Dwekl l2002). 
The use of exozodi in this context as has been referred broadly 
in the literature for Vega is misleading. 

The origin of hot dust population is unknown. One sce- 
nario that has been discussed frequently in the literature is 
that the hot dust arises from evaporating comets dynamically 
perturbed from the outer cold disks. In the case of Vega, a 
total mass of ~ 1O~ 9 M0 is requir ed to account for t he ob- 
served excess emission in near- IR dDefrere et al.ll20TTl) . The 
radiation blow-out time scale at 0. 1 AU is on the order of a 
year, suggesting a very high dust replenishing rate. This high 
replenishing rate implies that this phenomenon is unlikely to 
be in a static state. If not in steady-state, the hot dust may 
be created in transient dynamical events similar to the late 
heavy bombardment in our solar system, which was caused 
by a dynamical instability of the asteroid belt (Stro m et al.1 
2005). However, this kind of hot excess has also been found 
around stars that have no detectable infrared excess indica- 



tive o f the presence of a planetesimal population ( Abs il et ail 
2008), making this scenario unsatisfactory to explain the hot 

e xcess phenomeno n. 

Kobavash i et alj (12009) present an analytical model to form 
a narrow ring due to sublimation of dust grains drifting ra- 
dially inward due to P-R drag. This scenario could work in 
Vega and Fomalhaut if the dust drifting inward is refractory; 
i.e., the asteroid belt discussed in this paper being the source 
of the particles. This hypothesis only requires a P-R dom- 
inated disk; no dynamical perturbation is required. The re- 
plenishing rate implies that a total mass of > 0.4 M ffi in the 
Vega asteroid belt is required to retain such a rate over the 
age of the system (^400 Myr). However, the range of dis- 
tances and particle sizes from this drag-in component would 
re sult a much higher exc ess flux in the mid-IR (see Figure 5 
in lKobayashi et alJ201 ll) . Furthermore, one would expect that 
the amount of hot dust should be proportional to the amount 
of dust in the source region, i.e., the asteroid belt. The fact 
that the hot dust in both systems has a similar fractional lu- 
minosity, 5xl0" 4 (lAbsil et al.ll2006l l2009h while the warm 
dust in the Fomalhaut system is ^3 times more than that of 
the Vega system, argues against this non-stochastic, transport 
hypothesis. 

There may be an alternative direction for models of the 
hot dust component. The basic requirements are (1) a mech- 
anism that prevents the radiation-pressure blow-out of very 
small dust grains at very small astrocentric distances (~40 
stellar radii, i.e., the ^1500 K dust that is creating the K- 
band excess); (2) a grain composition that yields a feature- 
less, roughly Rayleigh- Jeans (or steeper than Rayleigh Jeans) 
spectrum between 2 and 10 /im (for consistency with the SED 
constraints in Section l3~TT i; and (3) a mechanism for the gen- 
eration grains of the appropriate composition. Although a 
detailed analysis is beyond the scope of this paper, these re- 
quirements evoke the behavior of nano dust grains (sizes of 
ten to a few tens of nm) that are charged by the solar wind 
and trapped by the solar magnetic field near the Sun. The un- 
derlying hypothesis is that sungrazing planetesimals deliver 
silicate-rich grains to the stellar neighborho od, but when the 
silica tes break down they yield metal oxides (Mann & Murad 
2005). Because of the high abundance of Mg and Fe in sil- 
icates, MgO and FeO are produced profusely in this process 
through reactions such as 



MgSiQ 3 ->MgO + SiO+-Q 2 



and 



Fe 2 Si0 4 



2FeO + SiO+-0 2 
2 



The resulting nano grains are very refractory, for example, 
MgO has a melting temperature of 3100 K and boiling tem- 
perature of 3873 K, nearly as high as the sublimation temper- 
ature of carbon, 3900 K; FeO boils at 3687 K. Consequently, 
MgO, FeO, and composite grains of the same elements can 
survive close to the star. We use the hot excess around Vega 
as an example to illustrate the resultant SED since the loca- 
tion for such a hot excess is best constrained among all the 
K band excess sources. We adopted the optic al constants for 
a mix ed iron/magnesium oxide, Mgo.eFeo 4O dHenning et al] 
1 19951) and computed their absorption and scattering efficien- 
cies using Mie theory. We computed a model SED (Figure 
ITOb for such a narrow ring between 0.18 and 0.2 AU from the 
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FIG. 10. — Model SED of the hot excess around Vega using nano oxide 
particles. Although the stellar spectrum shown on this plot (blue, thin, solid 
line) is the Kurucz atmospheric model (viewed pole-on), the actual heating 
spectrum is used the one viewed by the dust along the equator of Vega. The 
model spectrum (solid red line) and various power-law lines (as frequency v) 
are normalized to 2.2 fim to account for the 1.29% excess at K band. 



star, and composed of grains with radii of 5-20 nm (0.005- 
0.02 /im) in a size-distribu tion power index of -3.5 using the 
stellar model computed by Aufdenberg et al. (2006) (for dust 
viewed from the stellar equator). The dust temperatures for 
these gains in the region of 0.18-0.2 AU from Vega reach 
^2400 K. With the grain density of 4.8 g cm" 3 , a total mass 
of 2.3xlO _9 M ffi produces the 1.29% excess of the star at 2.2 
/im. Toward longer wavelengths to 10 /im, the output falls 
roughly as v 2 *, where Rayleigh Jeans falls as v 2 , that is the 
spectrum is somewhat steeper than Rayleigh- Jeans (see Fig- 
ure [Toll. There is a prominent feature at ~ 1 8 fim in our model 
SED due to the crystalline form of the material used to de- 
termine the optical constants, which might not represent the 
actual form of material resulted from sublimation of silicate- 
rich planetesimals. The mass in nano dust is equivalent to 
about 60 Halley Comet nuclei. Although the conversion of a 
comet nucleus to nano dust will not be fully efficient, it should 
be reasonably high since the process involves erosion of large 
grains and so most of the solid material (which is believed 
to constitute more than half the mass of a typical comet nu- 
cleus) may be converted into oxides. Thus, it appears that a 
plausible number of planetesimals (perhaps no more than a 
couple of hundred) would suffice to provide the nano grain 
population. Finally, we would like to emphyasize that the 
model SED shown in Figure [l0]is one example that satisfies 
the requirements (as listed previously) in our hypothesis. Our 
proposed scenario does not depend on the exact nature of the 
material as long as they meet these requirements. 

We now consider the magnetic trapping of these dust grains. 
Vega has a magnetic field strength of ~1.2 G, with about 
a quarter i n a dipolar component and the rest in higher or- 
der terms (Pet it et al.l 12.0101) . To first order, this is similar 
to that of the Sun, with a surface field of about 1 G, domi- 
nated by a dipole but with complex components due to ac- 
tivity. Nano grains in the vicinity of Vega will acquire elec- 
tric charge throug h the photoelectric effect. Thi s proce ss has 
been modeled by Pedersen & Go mez de Castrol (1201 ll) . who 
find that a 30 nm grain illuminated by an AO star will reach 
a level of about 200 e~, or charge per mass ratio (Q/m) ~2 
x 10~ 6 . Smaller grains will reach higher values of Q/m. As- 
suming a field similar to that of the Sun, these values are in 
the range where trapping occurs (ICzecho wski & Mann 2010). 
In summary, the hypothesis that the 2 /im excess is emitted by 



nano grains trapped in the magnetic field of Vega does sat- 
isfy the three conditions listed above. A better understanding 
of the behavior of the magnetic field of the star is required 
to improve our understanding of whether it is indeed the pro- 
cess that accounts for this emission. Furthermore, this mag- 
netically trapped, nano dust model can apply to other stars 
that show /T-band hot excesses through interferometry as long 
as these stars possess some magnetic field and can charge 
these nano particles through photoelectric effect (for early- 
type star) or stellar wind (for late-type stars). 

6. CONCLUSIONS 

Nearby debris disks, including the debris disk twins: Vega 
and Fomalhaut, have been playing an important role in our 
understanding of the underlying planetary architectures (plan- 
ets, and minor bodies) since their first discovery by IRAS 
through infrared excesses. Much attention has been focused 
on their large, cold Kuiper-belt-analog rings because they 
contain the majority of the left-over planetesimals and fine 
debris that covers a large surface area, making them readily 
detectable through infrared and submillimeter observations. 
The Fomalhaut debris system possesses an unresolved warm 
exce ss first suggested by t he Spitzer resolved image at 24 
/im (IStapelfeldt et al.ll2004T) . Using the Spitzer IRS spectrum 
centered at this warm excess along with the photometry mea- 
sured at the star position from the Herschel 70 /im and ALMA 
850 /im images, we corroborate that the warm (~170 K) ex- 
cess arises from thermal dust emission, an d is very unl i kely t o 
originate from stellar wind as suggested bv lAcke et al.l(l2012l) . 
Through resolved images at multiple wavelengths and mid- 
infrared spectrum, we identified for the first time the presence 
of a warm (~170 K) unresolved component in the Vega de- 
bris disk, which is clearly separated from the cold belt. Sim- 
ilar to the one in the Fomalhaut system, we suggest that this 
warm component also arises from thermal dust emission from 
a planetesimal belt located near the water-frost line, analo- 
gous to the asteroid belt in our solar system. The warm belts 
in both systems share many similar characteristics. No ex- 
tension at both MIPS 24 and PACS 70 /im (angular radius 
<3") and the observed dust temperature (~170 K) place this 
warm belt at ^2" from Vega and ~1."5 from Fomalhaut for 
blackbody-like grains. Furthermore, the stringent constraints 
on the excess levels from 3 to 6 /im around both stars indi- 
cate that this warm belt is not spatially associated with the 2 
/im (/T-band) excesses inferred from the interferometric ob- 
servations. 

Including the debris disks around e Eri and HR 8799, the 
four nearby debris systems all have spatially separated warm 
and cold belts where the observed dust temperatures for the 
warm belt are at ^150-170 K, whereas the cold belt tem- 
peratures are ^45-50 K despite the stark differences in stel- 
lar luminosity. Similar bi-modal dust temperature structure 
is also found in other debris disk systems. Systems where 
the location of the cold belts is known through resolved im- 
ages either in scattered light or thermal infrared/sub millimeter 
wavelengths and the association of a separate warm belt is in- 
ferred from spectrophotometric mid-infrared excesses include 
the following: HD 10647, HD 15115, HD 107146, and HD 
139664 (although not all of them have the warm-belt tem- 
perature at ^150 K). Bi-modal temperature distributions in- 
ferred from SEP ana lysis for many other unresolved sources 
(Morales eTaD|201JJ) argue that the debris disk structures are, 
somehow, determined by temperature-sensitive processes. It 
is well known that the location of the water-frost line in pro- 
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toplanetary disks plays an important role in the formation of 
planetesimals. An asteroid belt located near the frost line for 
these warm belts seems to imply that they are inherited from 
the early evolution of planetary systems. 

The /f-band excesses have been found around roughly a 
dozen of nearby main-seq uence stars using ground-based 
interferometric techniqu es (Absil et al. 200 6t Idi Folco et all 
120071; lAbsil etaill2008l 120091; lAkeson et al.ll2009t) . In some 
cases, sources like low-mass companions, ionized stellar wind 
or stellar scattered light have been ruled out except for hot 
(^1500 K) dust emission being responsible f or the excess 
emissi on. Recent A^-body simulations done by iBonsor et al.l 
d2012l) demonstrate the difficulty to scatter enough small bod- 
ies inward by a chain of planets inside a cold outer belt to sus- 
tain the observed level of small grains in Vega and r/ Crovi. 
In the case of Vega, the positive detections and constraints 
from non-detections using different interferometric facilities 
and wavelengths affirm the presence of very small grains in 
the very close vicinity (<0.2 AU) around the star. The prop- 
erties of these hot dust are very different from what we know 
about the particles in our own zodiacal cloud; instead, they 
are more likely to arise from the phenomenon observed near 
the Sun where nano particles trapped in the magnetic field 
of the star. In our proposed scenario, nano-size metal oxides 
originate from the sublimation of silicate -rich planetesimals, 
and are charged either via the photoelectric effect or the stel- 
lar wind, and then magnetically trapped in close proximity to 
the star. The replenishing rate of these tiny particles can be 
very low once they got trapped. Thus, our scenario does not 
require a massive reservoir of left-over planetesimals; there- 
fore, it presumably works for stars that show hot 2 /im excess 
but without detectable cold dust. 

Although the warm belts in the four systems discussed in 
this paper are not directly resolved, unlike the cold belts, the 
orbital ratios between the outer cold and inner warm belts are 
roughly >10 based on the observed dust temperatures. The 



large gap between the two planetesimal belts requires a mech- 
anism to maintain it mostly free of dust. The most plausible 
explanation is the existence of planetary bodies in the gaps, 
analogous to the HR 8799 and our own solar system. From 
simple chaotic zone calculations and the mass limit on Foma- 
lhaut b, we argue that multiple low-mass (<1 Mj) planets are 
required to maintain such a large gap. A similar argument is 
also supported for the Vega system, suggesting that the widely 
separated, two-belt debris systems are signposts for the pres- 
ence of multiple lo w-mass planets. O ur results are in line with 
the recent result bv lWvatt et al.l (1201 2|) where a positive corre- 
lation is identified between the debris disk detection rate and 
the presence of exoplanets with mass less than Saturn around 
60 nea rby G-type stars, and echo the recent R V and Kepler 
results dMavor et al.ll201 UlBatalha et al.ll2012l) that low-mass 
planets are more common than the massive ones, and that 
multiple planet systems are rather common among exoplan- 
etary systems. 
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